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Abstract 

The myxozoan Tetracapsu bides bryosalmonae is the causative agent of Proliferative Kidney Disease (PKD) targeting 
primarily the kidney of infected fish where it causes a chronic lymphoid immunopathology. Although known to be 
associated with suppression of some cellular aspects of innate immunity and a prominent lymphocytic hyperplasia, 
there remains a considerable knowledge gap in our understanding of the underlying immune mechanisms driving 
PKD pathogenesis. To provide further insights, the expression profiles of a panel of innate / inflammatory and 
adaptive immune molecules were examined in rainbow trout Oncorhynchus mykiss following a natural exposure to 
the parasite. Relative to controls, fish with early to advanced stages of kidney pathology exhibited up-regulation of 
the inflammatory cytokines interleukin (IL)-6 and IL-11, although remaining refractory towards genes indicative of 
macrophage activity. Antimicrobial peptides (AMPs) and anti-inflammatory markers, including cathelicidin (CATH) 
and IL-10 were markedly up-regulated during clinical disease. Up-regulation of adaptive immune molecules, 
including cell markers and antibody genes reflect the lymphocytic dominance of this disease and the likely 
importance of lymphocyte subsets in PKD pathogenesis. Up-regulation of T helper (T H ) cell-like response genes and 
transcription factors implies that T. bryosalmonae may elicit a complex interplay between T H cell subsets. This work, 
for the first time in the study of fish-myxozoan interactions, suggests that PKD pathogenesis is shaped by an anti- 
inflammatory phenotype, a profound B cell / antibody response and dysregulated T H cell-like activities. A better 
understanding of the functional roles of fish immune cells and molecules in PKD pathogenesis may facilitate future 
development of control measures against this disease. 



Introduction 

Proliferative Kidney Disease of salmonid fish is a slow 
progressive disease of major economic importance to 
aquaculture in the UK, Mainland Europe, and the USA 
[1,2]. Caused by the myxozoan parasite Tetracapsuloides 
bryosalmonae, PKD pathogenesis is a process closely 
linked to increasing water temperatures with climate 
change processes believed to impact on the latitude and 
altitude of the disease range [2]. Parasite spores, released 
from infected freshwater bryozoans, invade primarily via 
the gills and skin following the recognition of host nu- 
cleosides in mucous. Subsequently, parasites migrate 
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through the vascular system to organs including the kid- 
ney and spleen with the former being the main target 
organ for further development [3,4]. Extrasporogonic 
proliferation in the kidney interstitium provokes a 
chronic immunopathology characterized by a lympho- 
cytic hyperplasia, formation of granulomatous lesions, 
renal atrophy, and hyper secretion of immunoglobulins 
[1,2]. These stages are eventually eliminated in surviving 
fish and normal kidney function restored. Treatments 
with chemicals such as malachite green and fumigillin 
are effective against PKD, but these molecules are not 
licensable owing to their toxicity in humans [5,6]. Thus, 
in the absence of effective control measures immune 
therapy may be a way forward. However, this will neces- 
sitate an in-depth understanding of the immune 
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mechanisms underlying the kidney immunopathology 
and protective immune responses to PKD. This is par- 
ticularly poignant given that fish surviving PKD are 
immune to re-infection, providing the necessary impetus 
for vaccine development [1]. 

Previous studies examining host cellular responses 
during PKD have described a suppression of innate im- 
mune responses, including reduced phagocytic and re- 
spiratory burst activity of kidney phagocytes, reduced 
responsiveness to vaccination, and a dominance of im- 
munoglobulin (Ig) M-negative lymphocytes that undergo 
in situ proliferation [7]. Until recently, the availability of 
fish immune genes for gene expression studies was 
mainly limited to innate immune processes such 
as pattern recognition, antimicrobial protein, and com- 
plement activities, and cytokines involved in pro- 
inflammatory responses (e.g. IL-1(3 and tumour necrosis 
factor (TNF)-a) [8]. With the advent of fish genome and 
EST databases, there has been a tremendous increase in 
the sequencing and characterization of fish immune 
genes, particularly those homologous to cellular markers 
and response genes associated with CD4 + and cytotoxic 
T cell activity in mammals, providing tantalizing insights 
into fish T cell biology [8,9]. Likewise, the discovery of 
the fish-specific mucosal immunoglobulin class, IgT, has 
uncovered a subset of IgT-specific B cells, known to 
respond specifically to the intestinal myxozoan parasite, 
Cemtomyxa shasta and to the protozoan Ichthyophthi- 
rius multifiliis following gill invasion [10-12]. The chal- 
lenge now facing fish immunology is in the continued 
development of recombinant proteins and antibody 
markers, facilitating the functional characterization of 
fish immune responses and immune cell subsets in- 
volved in disease pathogenesis and protective responses. 
Given the current paucity of functional tools, expression 
profiling of fish immune genes represents a powerful 
means of providing insights into the immune mecha- 
nisms underlying disease pathogenesis and immune pro- 
tective responses, pinpointing potential areas for future 
immunological intervention [13-15]. 

Our previous studies have indicated that both T H i and 
T H2 -like immune processes are involved in PKD patho- 
genesis, whilst lacking the classical signs of a pro- 
inflammatory response in fish exhibiting different stages 
of clinical disease [16-18]. Given the dominance of 
proliferating lymphocytes/antibody levels during PKD 
pathogenesis, the known suppression of phagocyte activ- 
ity, the involvement of macrophages in the resolving 
stages of PKD and in response to other myxozoan para- 
sites, this study was undertaken to further examine 
genes indicative of innate/inflammatory and adaptive 
immune responses. Based on the current availability and 
functional characterization of rainbow trout immune 
genes, we targeted genes that encompass innate/ 



inflammatory/adaptive immune ligands, receptors, cell 
surface markers, antimicrobial peptides, and master 
transcription factors driving specific T H cell responses in 
higher vertebrates. Thus, our study provides a balance 
between potential innate and adaptive immune mecha- 
nisms shaping PKD pathogenesis. In line with chronic 
immune pathologies mediated by mammalian extracellu- 
lar parasites [19,20], our results suggest that PKD is 
characterized by a prevailing anti-inflammatory pheno- 
type, over-expression of immunoglobulin isotypes, and a 
dysregulated T H -like response. 

Materials and methods 

Fish sampling and monitoring 

Posterior kidney tissue samples were obtained from 
rainbow trout Oncorhynchus mykiss provided by a 
commercial trout farm in Southern England during a 
natural outbreak of PKD. The disease, at this site, is 
characterised by recurring annual epizootics with para- 
site detection in fish kidney smears occurring from late 
May to mid June depending on river temperature pro- 
files, flow rates, and other environmental factors [2,21]. 
The bryozoan population(s) harbouring T. bryosalmonae 
at this site is not known. However, based on the dynam- 
ics of parasite release from known bryozoan populations, 
it is estimated that parasites are present in the water 
course feeding the farm from mid March [2,21]. Two 
groups of fish from the same egg source {ca. 50-100 g 
each) were sampled for this study. Firstly, a parasite-na 
ive uninfected group and secondly, a parasite-naive 
group exposed to parasite-infected water from early 
April. Clinical signs of the disease were first seen in 
parasite-naive exposed fish early June. Sampling of both 
groups was undertaken late July at a water temperature 
of 15-16 °C when naive parasite-exposed fish exhibited 
kidney pathology ranging from early to advanced clinical 
stages (kidney swelling grades 1 to 3), as determined 
using the kidney swelling index system devised by 
Clifton-Hadley et al. [22]. All control fish had a kidney 
swelling grade of 0. In all fish sampled, approximately 
100 mg of kidney tissue was removed immediately below 
the dorsal fin, the area of the kidney associated with the 
onset of clinical disease. Tissue samples were placed 
into 1 mL of RNA-later (Sigma, ST. Louis, USA), kept at 
4 °C for 24 h and stored at -80 °C prior to RNA extrac- 
tion and PCR analysis. The presence of T. bryosalmonae 
kidney stages in parasite-exposed fish was confirmed by 
histological examination of posterior kidney smears and 
by qPCR (see below). Checks for other pathogen infesta- 
tions were conducted throughout the PKD season. With 
respect to fish used in this study, kidney swabs were 
streaked onto Tryptic Soy Agar (TSA) plates (Becton- 
Dickinson, Oxford, England) in order to check for 
the presence of common bacterial pathogens (e.g. 
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Aeromonas salmonicida). Streaked plates were incubated 
for 48 h at 20-22 "C prior to examination for bacterial 
growth. 

RNA extraction and PCR analysis 

Total RNA was extracted using TRI-reagent (Sigma) 
according to the manufacturer's instructions. Purified 
RNA was quantified using a Nanodrop spectrophotom- 
eter (NanoDrop Technologies, Wilmington, USA) and 
reverse transcribed into cDNA (20 ug per sample) using 
Bioscript (Bioline, London, UK) in 30 uL reactions. 
cDNAs were diluted to 500 uL with TE buffer (pH 8.0) 
and stored at -20 "C. For qPCR analysis, cDNAs 
representing uninfected fish (n = 7) and fish exhibiting; 
early (grade 1; n = 6), moderate (grade 1-2; n = 9), and 
advanced (grade 2; n = 10 and grade 3; n = 10) stages of 
clinical disease were examined. Sixty primer sets were 
used encoding putative cellular markers and immune re- 
sponse genes, as well as primers for the reference gene 
elongation factor- la (EF-la). A full list of all primers 
and associated information is provided in Additional file 1. 
EF-la has been repeatedly demonstrated to be highly 
consistent as a reference gene in the immune gene ex- 
pression profiling of fish host-pathogen interactions 
[13-15,23]. The relative parasite prevalence in all tissue 
samples was assessed by qPCR using T. bryosalmonae- 
specific primers for the detection of the house-keeping 
genes; 18S rDNA [EMBL: U70623] and 60S ribosomal 
protein L18: RPL18 [EMBL: FR852769] [24]. 18S rDNA 
primers were modified relative to those in previous PKD 
studies [25] and were used to detect the parasite in 
gDNA samples. Although T. bryosalmonae RPL18 is 
highly homologous to the rainbow trout homologue 
(53% amino acid identity), T. bryosalmonae-speciiic 
RPL18 primers were designed within the open reading 
frame aided by the considerable difference in codon 
usage between host and parasite (typically, 45-55% GC 
and 25-40% GC respectively). Hence, detection of T. 
bryosalmonae RPL18 transcripts provides a sensitive 
measure of only viable parasites in each sample, whereas 
18S rDNA detects both live and dead parasite material. 
With the exception of 18S rDNA, all primer pairs were 
designed and tested with a set of cDNA and DNA sam- 
ples to ensure that products could only be amplified 
from cDNA and not from genomic DNA under the con- 
ditions used. Genomic DNA was extracted from the 
same TRI-reagent tissue homogenates used for RNA 
extraction, as described previously [13]. SYBR green 
(Invitrogen, Paisley, UK) based RT-qPCR using Immolase 
DNA Polymerase (Bioline) was performed using a Light 
Cycler" 480 SW 1.5 system (Roche, Mannheim, Germany) 
as described previously [26]. For parasite DNA detection, 
primer efficiency was determined using serial dilutions of 
reference (internal PCR control) DNA and used for 



quantification of the DNA concentration. To normalize 
the level of T. bryosalmonae 18S rDNA to the input 
genomic DNA, additional qPCR was undertaken using 
primers to the trout macrophage colony stimulating 
factor (MCSF) gene, as described previously [13]. For 
cDNA detection, the primer efficiency and concentra- 
tion of each gene transcript was quantified using data 
generated from serially diluted reference DNA amplified 
in each PCR run, as described previously [13]. Since 
RPL18 mRNA detection is more sensitive than DNA de- 
tection and should more accurately reflect the viable para- 
site prevalence in individual fish, this was used for analysis 
of the immune gene expression data in addition to the 
kidney swelling grade assessment. To determine the RT- 
qPCR detection limit in terms of RPL18 transcript number, 
a pooled T. bryosalmonae positive sample was obtained 
from grade 2 cDNAs and serially diluted. A diluted RPL18 
reference was included to enable relative quantification. 
The expression of trout immune genes was initially nor- 
malized to the expression of EF-la and subsequendy 
expressed as fold change relative to the expression level in 
parasite-naive uninfected fish. Likewise, parasite RPL18 
cDNA levels were normalized to that of trout EF-la for 
each sample. 

Data analysis 

Relative immune gene expression levels were anchored 
to the lowest value of each data set and Log2 
transformed prior to statistical analysis, as described pre- 
viously [27]. Correlations between T. bryosalmonae 
prevalence (RPL18 mRNA detection), kidney swelling 
grade, and immune gene expression were assessed by 
calculating the Pearson product-moment correlation 
coefficient (r) and considered significant at P < 0.05 
(2-tailed). The significance of the average fold change 
between uninfected and infected groups was analysed by 
one-way analysis of variance (ANOVA) and the LSD 
post hoc test for comparison of means with differences 
considered significant at P<0.05 (2-tailed). All statistical 
analyses were performed using SPSS 8 Statistics package 
v 20.0 (IBM Corporation, Somer, New York, USA) and 
graphically represented using GraphPad Prism version 
5.04 (GraphPad Software Inc., La Jolla, USA). 

Results 

Pathogen detection 

Trout exposed to parasite-infected water were found to 
exhibit a range of clinical pathology ranging from grade 
1 at low levels to grade 3 in cases of severe/advanced 
kidney pathology. 18S rDNA levels steadily increased in 
kidney tissue samples from grade 1 to grade 3 fish, 
reaching maximal levels in grade 3 fish whilst remaining 
undetectable in uninfected fish (Figure 1A). RPL18 ex- 
pression, also undetectable in uninfected fish, plateaued 
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Figure 1 Detection of 7". bryosalmonae in rainbow trout kidney during clinical PKD. (A) T, bryosalmonae 18S rDNA relative to trout MCSF (DNA 
detection); (B) 1. bryosalmonae RPL18 relative to trout EF-1a (mRNA detection). Average cCq value in grade 1 fish is provided above the bar; (C) Correlation 
between 18S rDNA and RPL18 parasite detection markers is shown (r = 0.94; P< 0.0001); (D) Linear relationship of RT-qPCR detection and molar 
concentration of T. bryosalmonae RPL18. A serial diluted RPL18 internal reference was included in the same PCR run to enable relative quantification. 



at grades 1-2 and 2 decreasing by ca. 32% at grade 3 
relative to grade 2 (Figure IB). Nevertheless, parasite 
detection with 18S rDNA correlated highly with RPL18 
(r = 0.936), with both significantly correlating with swell- 
ing grade {r = 0.755 and 0.740 respectively) (Figure 1C). 
A linear relationship was apparent between RT-qPCR 
detection and RPL18 molar concentration with a detec- 
tion limit of 0.8 x 10' 7 M (ca. 20 transcript copies in a 
20 uL PCR reaction) (Figure ID). 

Importantly, agar plates prepared from kidney swabs 
from uninfected and T. bryosalmonae infected fish did not 
reveal the presence of bacterial pathogens, implying that 
the immune gene expression profiles observed in this 
study are influenced by the presence of T, bryosalmonae. 

Gene expression 

A prerequisite in assessing the biological significance of 
altered gene expression profiles is to consider the ex- 
pression levels indicated by average control Cq (cCq) 
values in uninfected fish. For example, cyclooxygenase 



(COX)-2B and inducible nitric oxide synthase (iNOS) 
transcripts in control fish are very rare (cCq: 32.7 and 
33.3 respectively) with a ca. 4- to 5-fold increase unlikely 
to be biologically significant, whereas a 3-fold increase 
in the highly abundant transforming growth factor 
(TGF)-|3la transcript (cCq: 22.9) is likely to be of greater 
biological significance. Therefore, in addition to fold 
changes in gene expression in grade 1 to 3 fish, the aver- 
age cCq values for the genes studied (i.e. in grade 0 fish) 
are presented in Table 1 and in Additional file 2. 

Overall, 36 gene transcripts correlated positively with 
parasite RPL18 transcript levels, with 9 genes (IL-6, IL- 
10A, IL-21, IgT (secretory), interferon (IFN) y, suppres- 
sor of cytokine signalling (SOCS)-3, CATH-1, CATH-2 
& hepcidin-1) strongly correlated (r>0.7) (Figures 2, 3, 
4, 5, 6, 7 and 8. See Additional file 2 for IL-18, IgT 
(membrane), type I IFN-A, IL-10B, SOCS-1, and IL-15 
data). Sixteen genes correlated positively with both para- 
site RPL18 and kidney swelling, including all genes 
strongly correlating with RPL18, together with IL-10B, 
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Table 1 Summary of immune gene expression in T. bryosalmonae infected rainbow trout kidney tissue 



Gene 


v_i 1 uuc \j 

average cCq 






Kidney swelling grade 




Grade 1 


Grade 1-2 


Grade 2 


Grade 3 


Pearson correlation 


T. bryosalmonae 18S rDNA 


Not detectable 


2105.69** 


4266.94** 


6652.34** 


8034.26** 


.755** 


T. bryosalmonae RPL18 




8013.47** 


37233.66** 


39385.81** 


26654.46** 


.740** 


IL-1B-1 


27.16 


3.71* 


2.24 


2.10 


2.31 


-.116 


TNF-al 


24.52 


1.12 


0.95 


0.94 


0.84 


-.155 


TNF-a2 


26.60 


0.97 


2.19* 


1.52 


1.64 


.106 


IL-6 


34.12 


7.44** 


31.53** 


23.59** 


1 5.38** 


.510** 


IL-11 


29.62 


31.74** 


25.04** 


14.13** 


11.54** 


.399** 


M17 


26.05 


5.10** 


6.28** 


5.48** 


2.60 


.179 


COX-2A 


24.69 


2.15 


1.68 


1 .75 


1.07 


-.161 


COX-2B 


32.73 


3.59 


3.94* 


3.10* 


2.03 


.164 


MCSF-1 


22.92 


1.18 


0.77 


0.49* 


0.23** 


-.658** 


MCSF-2 


21.78 


1.81 


3.82* 


2.77 


1.41 


-.076 


MCSF-R1 


27.65 


1.05 


2.48 


1.01 


0.57* 


-.375* 


MCSF-R2 


21.77 


2.15 


2.70 


2.65 


1.22 


-.001 


iNOS 


33.31 


2.34 


5.54* 


2.43 


1.28 


-.035 


Arginase-1 


15.20 


1.05 


0.50* 


0.45* 


0.23** 


-.612** 


Cathelicidin-1 


28.63 


52.60** 


136.21** 


72.80** 


57.81** 


.530** 


Cathelicidin-2 


29.58 


9.00** 


56.40** 


94.90** 


19.10** 


.547** 


Hepcidin-1 


25.76 


5.36* 


10.85* 


22.37** 


6.42** 


.360* 


LEAP-2A 


34.31 


10.92** 


12.87** 


7.76** 


6.84** 


.281 


CD8a 


25.31 


5.37** 


4.01** 


5.46** 


3.91* 


.238 


CD8B 


21.41 


2.99* 


3.05* 


3.34* 


2.33 


.128 


IL-2RB 


24.63 


3.79** 


5.19** 


1.28 


0.93 


-.168 


IgD H-secretory 


21.39 


0.98 


1.02 


0.89 


0.92 


-.095 


IgM H-secretory 


14.14 


5.39* 


3.38 


8.80** 


4.80** 


.399** 


IgT H-secretory 


21.31 


8.19* 


24.68** 


51.03** 


68.62** 


.732** 


CD4 


22.47 


3.66** 


4.28** 


4.72** 


3.56* 


.298 


T-bet 


24.94 


7.00** 


6.86** 


7.06** 


2.80** 


.203 


IL-2 


28.68 


3.68* 


4.16** 


2.16 


2.54 


.048 


IFNy 


29.25 


5.06** 


15.40** 


10.03** 


8.49** 


.500** 


GATA3 


25.00 


4.64** 


5.77** 


3.63* 


2.33 


.115 


IL-4/13A 


25.77 


1.67 


2.71 


1.70 


0.80 


-.124 


IL-21 


33.78 


20.28** 


42.62** 


31.58** 


11.91** 


.506** 


IL-22 


33.65 


10.03* 


6.67** 


5.49* 


4.48 


.159 


IL-17A/F2a 


28.27 


10.84** 


11.23** 


2.64 


3.24 


.061 


IL-17C1 


27.92 


6.56** 


2.74* 


1.43 


0.84 


-.259 


IL-17C2 


30.90 


2.08 


1.36 


1.06 


0.70 


-.349* 


IL-17D 


31.72 


4.95** 


5.73** 


3.78** 


3.56* 


.259 


FOXP3A 


26.19 


2.96 


5.15* 


4.75 


2.52 


.042 


FOXP3B 


27.03 


4.02** 


3.87** 


3.13** 


2.11 


.109 


TGF-BIa 


22.87 


3.34* 


3.00** 


2.02 


1.59 


-.015 
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Table 1 Summary of immune gene expression in T. bryosalmonae infected rainbow trout kidney tissue (Continued) 



nIL-IF 


24.76 


3.14* 


3.32** 


2.81* 


1.97 


.136 


IL-10A 


26.26 


14.68** 


18.21** 


25.24** 


9.51** 


.414** 


SOCS-3 


22.74 


3.44* 


8.38** 


10.24** 


5.90** 


.472** 



The average base-line Cq (cCq) in control (grade 0) fish is provided for each gene studied with the exception of parasite genes, which were undetectable in 
control fish. Immune gene expression at each stage of clinical PKD (grade 1 to 3} was initially normalized to trout EF-1a and subsequently expressed as fold 
change relative to expression levels in control fish. Pearson correlation r coefficients are given relative to kidney swelling grade. Significant differences (2-tailed) 
between control and infection groups are shown in bold. * P < 0.05, ** P < 0.01 . 




0 I ! . . ■ ■ 1 1 1 ■ » 0 I 0 

0 5 10 15 20 25 0 5 10 15 20 25 



T. bryosalmonae RPL18 prevalence T. bryosalmonae RPL18 prevalence 
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♦ 

0* ■ ■ 1 0 I I ■ ■ . , ■ ■ 
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T. bryosalmonae RPL18 prevalence T. bryosalmonae RPL18 prevalence 

Figure 2 Correlation of T. bryosalmonae prevalence and inflammatory genes in rainbow trout kidney. (A) IL-1 S-1; (B) TNF-a1; (C) TNF-a2, 
(D) IL-6; (E) IL-1 1; (F) M1 7. Pearson correlation r coefficients are given relative to T. bryosalmonae RPL18 and kidney swelling (in parentheses). 
Significant correlations are shown in bold. *P < 0.05; **P < 0.01 (2-tailed). 
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Figure 3 Correlation of T. bryosalmonae prevalence and macrophage marker/response genes in rainbow trout kidney. (A) COX-2A; 
(B) COX-2B; (C) MCSF-1; (D) MCSF-2; (E) MCSF-R1; (F) MCSF-R2. Pearson correlation r coefficients are given relative to T. bryosalmonae RPL18 and 
kidney swelling (in parentheses). Significant correlations are shown in bold. *P< 0.05; **P< 0.01 (2-tailed). 



IL-11, MCSF-1, arginase-1, IgM (secretory), IgT (mem- 
brane) and SOCS-1. Conversely, four genes correlated 
with the kidney swelling grade and not parasite RPL18, 
namely; MCSF-R1, RAR-related orphan receptor (ROR)y, 
IL-17C-2 and SOCS-2 (Figures 3 and 7. See Additional 
file 2 for RORy and SOCS-2 data). In each case a weak 
negative correlation was evident. Below, we focus on those 
genes correlating with parasite RPL18 expression, 55.6% 
(20 genes) of which would not have been detected by 



examining potential correlations between immune gene 
expression and kidney swelling grade alone. 

Expression of genes encoding the classical pro- 
inflammatory cytokines, IL-1|3 and TNF-oc did not cor- 
relate with parasite prevalence or kidney swelling grade 
regardless of the paralogues studied (i.e. IL-1|3-1 or IL- 
ip-3; TNF-ocl or TNF-a2). (Figure 2. See Additional file 2 
for IL-ip-3 data). In contrast, the inflammatory cyto- 
kines IL-6, IL-11, IL-18 and M17 all correlated with 
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Figure 4 Correlation of T. bryosalmonae prevalence and antimicrobial genes in rainbow trout kidney. (A) iNOS; (B) Arginase-1; 

(C) Cathelicidin-1; (D) Cathelicidin-2; (E) Hepcidin-1; (F) LEAP-2A. Pearson correlation r coefficients are given relative to T, bryosalmonae RPL18 
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parasite prevalence, with IL-6 markedly so (r = 0.85). IL- 
6 and IL-11 gene expression also correlated with kidney 
swelling grade (Figure 2). Another proinflammatory me- 
diator, COX-2, also showed a significant correlation with 
parasite prevalence (Figure 3), but only with respect to 
COX-2B. Both MCSF-1 and arginase-1 exhibited a sig- 
nificant negative correlation with both parasite preva- 
lence and kidney swelling grade (Figures 3 and 4). The 
expression of both isoforms of the MCSF receptor 
(MCSF-R1/2) did not correlate with parasite prevalence 



and only weakly correlated with kidney swelling grade in 
the case of MCSF-R1 (Figure 3). Similarly, the expression 
of a marker of mature antigen presenting cells (CD83) did 
not correlate with either infection read-out. Gene expres- 
sion of the antimicrobial peptides examined correlated 
with parasite prevalence, strongly in the case of CATH-1/2 
and hepcidin; furthermore, these genes also correlated 
significantly with kidney swelling grade (Figure 4). 

To examine the potential contribution of the adaptive 
immune system to PKD pathogenesis, the expression of 
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genes homologous to mammalian B and T cell markers 
and signature molecules of T H responses were investi- 
gated. Expression of secretory forms of IgM and IgT 
correlated positively with parasite prevalence and 
kidney swelling grade, with IgT correlating strongly 
(r= 0.747 and 0.732 respectively) (Figure 5). The 
membrane form of IgT also exhibited a positive, albeit 
weak, correlation with both infection read-outs, whilst 
membrane IgM was refractory. In contrast, membrane 
and secretory forms of IgD were refractory to both 



infection read-outs with no PKD-mediated gene 
modulation detected (Figure 5 Additional file 2). Ex- 
pression of the T cell markers CD4 and CD8 (sub- 
units a and (3) exhibited a positive correlation with 
parasite prevalence but not with kidney swelling 
grade. Strikingly, the IL-2 receptor chain, IL-2R|3 
although not correlating with either infection read- 
out was up-regulated at early clinical stages, dropping 
sharply at grades 2 and 3, as shown in Figures 5 and 
6, and Additional file 2. 
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Significant correlations are shown in bold. *P < 0.05; **P < 0.01 (2-tailed). 



In depth analysis of genes encoding T H i signature 
molecules revealed a putative T H i-like component to 
PKD pathogenesis. The master T H1 T cell-specific T-box 
transcription factor, T-bet, was positively correlated with 
parasite prevalence (Figure 6). The expression of IFN-y 
and IL-2 correlated with parasite prevalence, with IFN-y 
also correlating with kidney swelling grade (Figure 6). As 
shown in Additional file 2, another important growth 
factor for T cells, namely IL-15, was upregulated during 



early clinical stages, correlating positively with parasite 
prevalence. The master T H 2 trans-acting T cell-specific 
transcription factor, GATA3, exhibited a similar expres- 
sion profile to T-bet (Figure 6). However, expression of 
the trout IL-4/13 cytokine family member, IL-4/13A, 
was not affected by infection (Figure 6). In the context 
of T H17 cell-derived cytokines, we have analysed the 
transcriptional profiles of the four IL-17 family members 
found in trout, namely, IL-17A/F2a, IL-17C-1/C-2, and 
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IL-17D (Figure 7). IL-17A/F2a, IL-17C-1, and IL-17D 
exhibited only weak or no correlation with parasite 
prevalence and no correlation with kidney swelling 
(Figure 7). Strikingly, however, in a similar manner to 
IL-2R(3, IL-17A/F2a and IL-17C-1 were both upregulated 
at early clinical stages, dropping sharply at more 
advanced stages (Table 1). Expression of the trout 
homologue of the T H17 master transcription factor, 
RORy did not correlate with either infection read-out. In 
contrast, the T H17 cytokine, IL-21 strongly correlated 



with parasite prevalence (r =0.737) and significantly, 
although to a lesser extent, with kidney swelling grade 
whereas, IL-22 correlated only with parasite prevalence 
(Figure 7). 

A range of anti-inflammatory immune genes was also 
studied, including IL-10 paralogues, TGF-pia, the trout 
IL-1|3 antagonist nIL-lF, paralogues of the master T reg 
transcription factor, forkhead box (FOX) P3, and several 
members of the SOCS family (Figure 8. See Additional 
file 2 for SOCS-1, -2, -5b, -7 data). The expression of 
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both IL-10 paralogues strongly correlated with parasite 
prevalence (IL-10A: r = 0.735 and IL-10B: r = 0.679) and 
weakly with kidney swelling grade, whilst the expression of 
TGF-pia, nIL-IF, and FOXP3A/B exhibited weak or no 
correlation with both infection read-outs (Figure 8. See 
Additional file 2 for IL-10B data). Of the SOCS molecules 
examined, SOCS-1 and SOCS-3 correlated significantiy 
with parasite prevalence and kidney swelling grade 
(Figure 8. See Additional file 2 for SOCS-1 data). Lastly, 
four other genes of relevance to inflammatory processes 



were studied, namely CD9, type I IFN, vascular endothelial 
growth factor (VEGF), and ciliary neurotrophic factor 
(CNTF). As shown in Additional file 2, type I IFN 
exhibited a weak correlation with parasite prevalence, but 
no correlation with kidney swelling grade. All other genes 
were refractory towards both infection read-outs. 

Discussion 

Until recently, the paucity of available fish immune genes, 
particularly those indicative of T H cell-like activities, has 
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hampered gene expression studies investigating fish host- 
pathogen interactions. Here, for the first time in the study 
of host-myxozoan interactions, we present gene expres- 
sion data providing insights into the underlying innate and 
adaptive immune mechanisms shaping the pathogenesis 
of PKD in trout. Such knowledge may pinpoint immune 
molecules and pathways that could be targeted for im- 
munological intervention. 

In general, immune gene expression changes were 
more strongly correlated with parasite prevalence than 
kidney swelling. With respect to parasite prevalence, the 
apparent reduction in viable parasites at swelling grade 3 
from the parasite RPL18 data relative to 18S rDNA 
could signify the activation of mechanisms leading to 
pathogen clearance and recovery. Differences between 
parasite prevalence and kidney swelling grade are likely 
to be due to the fact that, for each fish, parasite preva- 
lence and immune gene expression data were collated 
from the same tissue sample, whereas kidney swelling 
data is based on a qualitative holistic assessment of 
disease progression. The modulation of nine genes cor- 
related strongly with parasite prevalence with IL-6, 
CATH-1, and CATH-2 exhibiting r values > 0.8. Recent 
studies describing functional associations between im- 
mune gene expression and parasite prevalence in wild 
host populations suggest that they may reflect host- 
parasite co-evolutionary processes [28]. Similarly, trout 
immune genes strongly correlating with T. bryosalmonae 
prevalence in the current study may exemplify host 
genes under selection. 

Innate immunity / inflammatory mediators 

The cytokines, IL-lp and TNF-oc, are crucial mediators 
of pro-inflammatory responses and in the activation of B 
and T cells [29]. Functionally assessed to have pro- 
inflammatory activity in fish, IL-1[3 and TNF-a are often 
co-expressed with other macrophage-derived inflamma- 
tory mediators such as IL-8, COX-2, and iNOS in parasitic 
and bacterial infections [13-15,30,31]. Our previous stud- 
ies revealed no or a very modest impact of PKD on the 
gene expression of pro-inflammatory mediators, which 
has been corroborated by the current study [16]. Similarly, 
other myxozoan fish parasites appear to elicit either weak 
transient up-regulation or, indeed, down-regulation of 
genes encoding pro-inflammatory mediators as part of a 
global down-regulation of innate and acute phase response 
genes [32,33]. 

In the current study, we extended the repertoire of 
inflammatory mediators to include IL-6, IL-10, IL-11, 
TGFpl, and the IL-1[3 antagonist, nIL-lF. Intriguingly, 
IL-6 and IL-11 transcripts were up-regulated over 
30-fold correlating positively with both infection read- 
outs. Unlike IL-1|3 and TNF-a, IL-6 exhibits both pro- 
and anti-inflammatory activities in the regulation of 



macrophages and lymphocytes [34]. Trout studies have 
also attributed pro- and anti-inflammatory activities to 
fish IL-6 in stimulating macrophage proliferation and 
induction of cytokine inhibitors (SOCS-1/-3), antimicro- 
bial peptides (CATH-2, hepcidin, & liver-expressed 
antimicrobial peptide-2 isoform A; LEAP-2A), whilst 
down-regulating TNF-a and IL-1[3 [35]. Furthermore, 
PKD elicited transcriptional up-regulation of the fish- 
specific IL-6 family member, M17 which correlated 
positively with parasite prevalence [36]. We previously 
described PKD-associated up-regulation of SOCS 1/3 
genes [18]. Thus, marked up-regulation of CATH-1/-2 
and LEAP-2A along with the unresponsiveness of TNF- 
a and IL-1[3 in the current study could be, in part, due 
to the activity of trout IL-6 family members. Fish IL-10 
has been shown to be anti-inflammatory, suppressing 
pro-inflammatory cytokine expression, phagocytosis, and 
respiratory burst activity [37]. Considering the domin- 
ance of IL-10 gene expression during PKD pathogenesis, 
this cytokine could be an important player in the 
PKD-mediated suppression of phagocytic activity [7,37]. 
Interleukin 11, also markedly influenced during PKD, is 
considered to be a potent anti-inflammatory cytokine in 
mammals [38]. Although not functionally characterised in 
fish to date, IL-11 could, nevertheless, also play an import- 
ant role in the suppression of pro-inflammatory mediators 
and innate immunity during PKD pathogenesis. 

The apparent lack of a pro-inflammatory response to 
T, bryosalmonae at the transcriptional level is particu- 
larly intriguing given that fish bacterial and other 
parasitic infections exhibit, to some extent, concomitant 
up-regulation of pro-inflammatory (TNF-a, IL-1[3, IL-8) 
and anti-inflammatory (IL-6, IL-10, IL-11, nIL-lF) cyto- 
kines [13-15]. This disparity could be attributed to the 
characteristic lymphocytic dominance of PKD pathogen- 
esis relative to other fish diseases, although macrophages 
are as prominent as lymphocytes in granulomatous 
lesions that develop during the resolving stages of PKD 
[1]. The putative lack of a macrophage-driven pro- 
inflammatory response to T. bryosalmonae is further 
reinforced by the refractoriness or, indeed, down- 
regulation of other macrophage marker genes such as 
arginase-1 and MCSF ligand/receptor paralogues [39]. 

A surprising finding in this study was the strong 
correlation between parasite prevalence and the promin- 
ent transcriptional up-regulation of the antimicrobial 
peptides, CATH-1, CATH-2, and hepcidin-1. Antimicro- 
bial peptide expression is known to be influenced by 
parasite infection in both vertebrate and invertebrate 
host-parasite interactions. Salmonid CATH-1 and/or 
CATH-2 genes are markedly up-regulated by both bac- 
terial and parasitic infection [14,23,35]. CATH-1/-2 pep- 
tides have been ascribed direct antibacterial and indirect 
cytokine modulating activities in fish, whilst a CATH-2 
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peptide, was able to slow oomycete sporulation [14,40]. 
These studies highlight the potential for fish cathelici- 
dins to protect against T. bryosalmonae directly or indir- 
ectly through the modulation of cytokine activity. As in 
fish, mammalian cathelicidins are influenced by bacteria 
and parasites mediating both pro- and anti-inflamma- 
tory responses. Although considered beneficial, they may 
also contribute to infection-mediated immunopathology 
[41]. Thus, cathelicidins may add to the apparent anti- 
inflammatory nature of PKD and to the associated patho- 
genesis. Hepcidin, an iron regulator that reduces the 
availability of iron for pathogen growth, is up-regulated by 
bacterial and parasitic infection in salmonids [14,23,42]. 
Although not functionally characterized in fish, an import- 
ant role for this molecule in the control of mammalian 
parasitic infections has been recently demonstrated [43]. 
As with cathelicidins, fish hepcidin is likely to have a 
broad spectrum of activity against fish pathogens. 

Adaptive immunity 

Suppression of innate immune/phagocytic activity and 
in situ proliferation of lymphocytes are two of the most 
significant cellular aspects of clinical PKD [7]. An over- 
shadowing anti-inflammatory phenotype may account 
for the aberrant lymphocytic response through suppres- 
sion of phagocyte activity. Hyperimmunoglobulaemia is 
also a major characteristic of PKD with immunoglobu- 
lins also playing a major role in other myxozoan infec- 
tions [1,10]. Here we report the presence of very high 
levels of transcripts encoding the pathogen-responsive 
immunoglobulin isotypes, IgM and IgT. All Ig primer 
sets were designed in the heavy chain constant region. 
We, therefore, acknowledge that we were unable to 
distinguish between viable and sterile Ig transcripts. 
However, we know from recent trout studies that the 
abundance of viable Ig transcripts during infection pro- 
cesses is ca. 70% [44]. Intriguingly, the membrane form 
of each Ig subtype was not or very weakly influenced 
whereas the secretory form was markedly up-regulated, 
which agrees with the IgM negative nature of the prolif- 
erating lymphocytic population in the kidneys of T. 
bryosalmonae-infected fish [7]. Our data corroborates 
with previous reports of hyperimmunoglobulinaemia in 
that Ig secretion appears to be a major feature of PKD 
pathogenesis [1]. The absence of up-regulated mem- 
brane Ig transcripts may suggest that B cells are not 
undergoing proliferation during the course of PKD and 
that the IgM negative proliferating lymphocyte popula- 
tion associated with PKD pathogenesis could be attrib- 
uted to the proliferation of T cell subsets. This would, 
however, need to be further substantiated using specific 
antibodies to trout IgT and T cell markers when avail- 
able. In support of this premise, we observed significant 
transcriptional up-regulation of CD4 and CD8 sugges- 



ting a potential role for T H and T c -like responses in 
PKD pathogenesis. Aberrant T cell activity, could to 
some extent account for the over-expression of antibody 
isotypes as seen in mammalian parasite-mediated chro- 
nic immunopathologies [20]. 

To investigate T cell-like activity in PKD, the ex- 
pression of trout genes homologous to mamma- 
lian molecules defining T H i, T H 2> T H i7, and T reg -like 
responses were studied [9]. Mammalian pathogen- 
mediated chronic immunopathologies have been linked 
to dysregulated T H responses involving a complex inter- 
play between stimulatory and suppressive immune sig- 
nals leading to pathology rather than parasite clearance 
[19,45]. IL-17A is one of the main drivers of parasite- 
mediated pathology, particularly in the formation of 
granulomatous tissue, a process regulated by IL-10, 
TGF-(3 and T reg cell activity. Reducing IL-17A in this 
context enhances protective mechanisms by reducing 
pathology/parasite prevalence and enhancing parasite- 
specific antibodies [19]. Furthermore, T reg responses are 
thought to be responsible for T cell anergy under such 
conditions resulting in reduced responsiveness to by- 
stander vaccines [20]. PKD also entails the formation of 
granulomatous lesions and reduced responsiveness to 
by-stander vaccines [7]. Thus, fish T H -like activities may 
play an integral role in shaping this disease. In line with 
previous fish studies, T-bet, the master T H i transcription 
factor and the T H i signature cytokines, IFN-y and IL-2 
(including receptor) were up-regulated during PKD 
[13,15,17]. Both cytokines are known to drive T H i-like 
responses in fish, which implies that in the context of 
adaptive immunity, PKD appears to elicit T H i-like activ- 
ities [46,47]. T H2 -like responses were also evident as 
judged by up-regulation of the T H 2 master transcription 
factor, GATA3, whereas the trout T H 2-like cytokine, IL4/ 
13A was refractory towards clinical PKD [48,49]. Other 
trout IL4/13 isoforms yet to be discovered, however, 
may play a more prominent role in clinical PKD. The 
T H i7 master transcription factor RORy (a splice variant 
of RORy) has not been described in fish to date, al- 
though isoforms of RORy have been sequenced and 
characterized in rainbow trout [50]. The expression of a 
trout RORy transcript did not correlate with either in- 
fection read-out. As with IL-4/13A, the existence of 
more PKD-responsive gene paralogues or other splice 
variants may also account for the apparent refractoriness 
of RORy and the modest changes of trout IL-2 [47,51]. 
The T reg master transcription factor, FOXP3 exists as 
two isoforms in trout both of which are up-regulated by 
the T cell mitogen, Phytohaemagglutinin (PHA) [48]. 
Here, both isoforms were transcriptionally up-regulated 
during clinical PKD, which may suggest a role of trout 
T reg -like cells in this disease, particularly as the cytokine 
mediators of FOXP3 expression, IL-10 and TGF-|3la, 
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were also up-regulated. Evidence for T H i7-like activities 
in fish has been provided in the sequencing and 
characterization of four IL-17 genes, all of which have 
been shown to respond to bacterial infection and other 
immune insults, with trout IL17A/F2a demonstrating 
pro-inflammatory activity [8,52]. As with IL-17, IL-21 
and IL-22 are primarily derived from T H i7 cells and are 
integral to the functioning and proliferation of the T H i7 
phenotype [53,54]. Initial differentiation of T H17 cells is 
driven by IL-6 and TGF-(3 and subsequently by IL-21 
and TGF-[3 [19]. Functional studies have confirmed trout 
IL-21 as a key regulator of B and T cell-like activities in- 
ducing IFN-y, IL-10 and IL-22, and maintaining CD8 
and IgM transcript levels in head kidney leucocyte cul- 
tures and increasing CD4, T-bet and GATA3 transcripts 
in the same cells [27]. Here, clinical PKD was associated 
with marked transcriptional up-regulation of IL-21, 
IL17A/F2a, and to a lesser extent, IL-22 and IL-17C-1, 
although only IL-21 correlated positively with parasite 
prevalence and kidney swelling. These data, in conjunc- 
tion with up-regulated IL-6 and TGF-[31a suggests that 
the lymphocytic character of PKD could involve a T H i7- 
like activity. 

In conclusion, our results imply that PKD pathogen- 
esis is shaped by a profound B cell/antibody response 
and potentially dysregulated T H cell activity. Our data is 
also suggestive of phagocyte-mediated pro-inflammatory 
processes being over-shadowed by a prevailing anti- 
inflammatory phenotype during PKD pathogenesis, al- 
though T H1 -like mechanisms do appear to be stimulated 
by PKD. Although many T cell marker and response 
gene homologues have been identified in fish, there is 
still a current knowledge gap in terms of the functional 
characterization of fish T cells. Specific antibodies to 
immune cell markers and neutralizing antibodies are 
required to aid future functional assessment of the rela- 
tive contribution of immune cell populations and im- 
mune regulatory molecules to the pathogenesis of PKD 
and other fish diseases. Such tools will greatly facilitate 
the identification of immune therapies to drive effector 
responses towards host protection rather than chronic 
immunopathogenesis. 

Additional files 



immune gene expression. Statistical data (P-values) representing gene 
expression (fold change) at each swelling grade and correlation analysis 
of gene expression relative to T. bryosalmonae RPL18 and kidney swelling 
grade are given. Significant differences (2-tailed) are shown in bold. 
*P<0.05, **P<0.01. 
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